p self energy at finite temperature and density in the real-time formalism 



Sabyasachi Ghosh and Sourav Sarkar 
September 14, 2011 

7— I ' 

o 

CnJ ■ Theoretical Physics Division, Variable Energy Cyclotron Centre, 

p jl 1/AF, Bidhannagar, Kolkata 700064, India 

Abstract 

The p meson self-energy in nuclear matter from baryonic loops is analysed in the real time formulation 
of field theory at finite temperature and density. The discontinuities across the branch cuts of the self- 
energy function are evaluated for an exhaustive set of resonances in the loops considering the fully relativistic 
thermal baryon propagator including anti-baryons. Numerical calculations show a significant broadening of 
the p spectral function coming from the Landau cut. Adding the contribution from mesonic loops, the full 
spectral function of the p in a thermal gas of mesons, baryons and antibaryons in equilibrium is evaluated at 
various values of temperature and baryonic chemical potential. 

1 Introduction 

> : 

■ It is well known that the rate of dilepton production from a thermal system is proportional to the two-point 
£ — . correlation function of vector currents. Hence the spectral function of vector meson, the p meson particular, plays 
such an important role in the analysis of the late stages of heavy ion collisions [HE]- The NA60 experiment at 
Q\ ' the CERN SPS measured dimuon pairs in In-In collisions in which an excess was observed over the contribution 
from hadronic decays at freeze-out in the mass region below the p peak [3] . This was attributed to the broadening 
of the p in hot and dense medium [2]. More recently, the PHENIX experiment reported a substantial excess of 
electron pairs in the same region of invariant mass [4]. This has been investigated by several groups but the 
yield in all these cases have remained insufficient to explain the data. Thus the issue of low mass lepton pair 
yield in heavy ion collisions is far from closed and is one of the key issues to be addressed in the forthcoming 
& \ Compressed Baryonic Matter (CBM) experiment to be performed at the FAIR facility in GSI [5]. 

A substantial volume of work has been devoted to the study of p meson properties in hot and dense medium. 
We do not attempt to review the existing literature but mention a few of them to put our work in perspective. 
We find that it is only for the tt — ir loop [6l [7] that one calculates the thermal loop directly. In the case of other 
loops typically involving one heavy and one light particle or both heavy particles one uses in general either 
the virial formula [U [10] or the Lindhard function \12\ [T3| [T4"l [T5] . Most of the calculations involving 
baryonic effects mentioned above were performed at zero temperature. Finite temperature effects on the p 
spectral function in dense matter have been evaluated by Rapp et al [12] in terms of resonant interactions of the 
p with surrounding mesons and baryons in addition to modifying the pion cloud. Eletsky [10] and collaborators 
have also evaluated the spectral function of vector mesons at finite temperature and density in terms of forward 
scattering amplitudes constructed using experimental inputs assuming resonance dominance at low energies and 
a Regge-type approach at higher energies. 

The sources modifying the free propagation of a particle find a unified description in terms of contributions 
from the branch cuts of the self energy function as shown by Weldon [16] . In addition to the unitary cut present 
already in vacuum, the thermal amplitude generates a new cut, the so called the Landau cut which provides 
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the effect of collisions with the surrounding particles in the medium. This formalism was applied to obtain the 
p self-energy in hot mesonic matter |17j by evaluating the one loop self-energies involving the tt, uj, hi and a\ 
mesons. A significant broadening of the spectral function was obtained without appreciable shift in the mass 
as expected from chiral interactions. 

In this work, we extend this analysis to the case of baryonic matter at finite temperature considering an 
exhaustive set of 4-star resonances in the baryonic loops making up the p self-energy. The framework of real 
time thermal field theory [18j 113 US EI] that we use, enables us to evaluate the imaginary part of the self-energy 
from the branch cuts for real and positive values of energy and momentum without having to resort to analytic 
continuation as in the imaginary time approach [22]. Here we work with the full relativistic baryon propagator 
in which baryons and anti-baryons manifestly appear on an equal footing. Thus the contributions from all the 
singularities in the self-energy function including the distant ones coming from the unitary cut of the loops 
involving heavy baryons are also included. These are usually not considered but can contribute appreciably to 
the real part of the p meson self-energy as shown [23] in the case of a NA loop. In addition we have used the 
covariant form of the momentum dependent vertex functions in the loop integrals in which additional terms |24j 
required to describe the coupling of off-shell spin 3/2 fields have been introduced. 

In the following section we define the correlation function of vector currents and its relation to the transverse 
propagator of the p. We also provide the various Lagrangian densities which will be used at the vertices of the 
loop graphs. Next, in section 3 we specify the kinematic decomposition of the thermal propagators. In section 
4 we evaluate the baryonic self-energy graphs as well as the discontinuities across the branch cuts. Section 5 
contains the results of the numerical evaluation of the real and imaginary parts as well as the p spectral function 
followed by a summary and discussions in section 6. In the appendix we provide the details of various factors 
appearing in the expression for self-energy for the different loops and provide some details of evaluation of the 
imaginary part in addition to a brief discussion on propagators and self-energies in the real time formalism. 



2 The two point function in the medium 

We begin our discussion with the two point function of vector currents in vacuum, 

J%(E,q) = i J d 3 xdre^(0|Ty;(x)^(0)|0) (1) 
where VHx) are the vector currents of two flavour QCD, given by 

V*(x) = q(xh^q(x), 1=( U d ) ( 2 ) 

t 1 being the Pauli matrices. In the real time formulation of thermal field theory, the in medium two point 
function assumes a 2 x 2 matrix structure [2D]. The thermal two point function is given by 

Tfi ab (E,q) = i J d 3 xdTe iq x {T c V*(x)Vj(0)) ab (3) 

where (O) denotes the ensemble average of an operator O, 

(O) = Tr{e~ pH 0)/Tre~ pH (4) 

and Tr indicating trace over a complete set of states. The superscripts a, b (= 1, 2) are thermal indices and T c 
denotes time ordering with respect to a contour in the plane of the complex time variable [21]. The two point 
function of vector currents can be related to the p meson propagator using the method of external fields [25] 
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where one introduces a classical vector field v^x) coupled to the vector current V p (x). The free propagator of the 
rho meson can be obtained by coupling the external field to the p meson field operator using the Lagrangian [26] 

p 

C P v = —d^v" • (<9 p„ - d u p p ) 
m p 

where F p = 154 MeV is obtained from the decay p° — > e + e~ . 

The transverse p meson propagator G a pv is then obtained from the relation Tj^ = K p G a pv where the factor 
K p = (F p q 2 /m p ) 2 comes from the coupling of the current with the p field [17]. The isospin structure is given 
by 5 lJ which we omit from now on. 

The free propagation of the p meson is modified by interactions in the medium which is populated by mesons 
and baryons. Here we consider one loop graphs shown in Fig. [T] consisting of the nucleon N and another baryon 
R denoted by the double lines. We have included all spin one-half and three-half 4— star resonances listed by 
the PDG [27] so that R stands for the iV*(1520), JV*(1650), iV*(1700), A(1230), A*(1620), A*(1720) as well as 
the iV(940) itself. Omitting isospin factors, the pN couplings with the resonances are described by the gauge 
invariant interactions [15] 

/ — -uu . .i. , l _ i t p _ 1 + 



C = — [iiRO-^PuvipN + h.c] J 



R 
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C = — [ip R Y7 PpuipN + h.c] J\ 



R 
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where p pu = d p p u — d u p p and = §[7^7^ — 7 V 7 M ]- The isospin part of the RNp interaction is given by 

\^Rabc(r-p) b d r N e cd 1 = 3/2 (6) 



V2 

where the indices a,b,c,d take values 1 and 2 and e 12 = — e 21 = 1. For the self-energy diagrams shown in Fig 
1, the isospin factor Ip comes out to be 2 for / = | and | for I = |. 

It is essential to point out that for the spin 3/2 resonances this coupling is not quite correct owing to the fact 
that the free Lagrangian for the Rarita-Schwinger field ip^ has a free parameter [28]. A symmetry is associated 
with a point transformation under which the free Lagrangian remains invariant up to a change in the value of 
the parameter [23]. The standard practice is to make a choice of the value of this parameter so that the spin-3/2 
propagator has a simple form. In order that the interaction also remains invariant under this transformation 
an additional term is added to it. Thus the Lagrangians involving spin-3/2 fields take the form 

£ = -L$ a R O a ^y^ N + h.c.] j£ = ^ (7) 

with O pa = g pa — j "f p 7c, the second term contributing only when the spin 3/2 field is off the mass shell. The 
value of the coupling strength / thus remains unaffected by this exercise. 
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Figure 1: One- loop Feynman diagrams for the two-point function contributing to the p self-energy in baryonic 
matter. The solid and double lines stand for nucleons and resonances respectively. 

3 Kinematics of the p propagator 

The complete propagator of the p is obtained from the Dyson equation [201 ES] 

where n{^.' ab denotes the thermal self-energy matrix and G^} ah {q) stands for the free thermal propagator. 

As described briefly in the appendix, one can get rid of the thermal indices by diagonalisation. In terms of 
the diagonal elements (denoted by bar) which are analytic functions, the Dyson equation for the p propagator 
reads 

GV(g) = G^{q) -G^(q)T^: t ( q )G au (q) , (9) 

where 



(q) = -9,u + ~ 2 J— • (10) 

F V q J q — m P + ie 

The one loop self energy with baryons is obtained from the two diagrams shown in Fig. Q] so that 

tig( q ) = tf"( q ) +n a \- g ). (11) 

On addition of the contribution from the meson loops the total p self-energy is given by 

ni: t (q) = n X B(q) + TiM(q) • (12) 

In the medium, the presence of the four velocity introduces an additional scalar variable u ■ q in addition to 

q 2 leading to two independent tensors P^ and in terms of which the propagator and self-energy can be 
written as 

Gfj, v = PpuGt + Q^i/Gi 

IV = P^Ut + Q^Ut (13) 



with 



D . QiaQv q 2 ~ ~ ~ i v , 2 

P/iv = -9/iv + -^2 =2 u v u vi u li = u ll - [u- q)q v ,lq ; 

Qiiv = —z^-u^Uy, q 2 = (u- q) 2 - q 2 . (14) 



Using (fP3"|h the Dyson equation @ can be solved to get, 

-1 

m 2 -U t {qy q 2 q 2 - m 2 - q 2 Ili(q) 



GM = , G t (q) = - 2 - ^ (15) 



4 



where 



1 ,— u Q 2 ^ s TV i ^ TT 

2 



n* = —(nj + ^n 00 ), n, = -^n 00 , n 00 = uwn^ . (16) 



The self-energy function can be obtained from the 11-component of the in- medium self-energy matrix using 
(see appendix) 



RelV = Rell^ 



Imn p = e(g )tanh(/3g /2)Imn^ (17) 



in terms of which the retarded self-energy is given by 

Re R^j, = Re 

ImlV = e(go)Imn^ . (18) 
We now proceed to evaluate the 11-component of the rho self-energy in the following section. 

4 The self energy and its analytic structure 

Let us begin by writing the expression for the p self-energy in vacuum corresponding to the first diagram in 
Fig. 1. For spin 1/2 resonances in the loop, this is given by 



IF-(g) = U F (t^j 2 J ^Tr[T»S(p, m N )T»S(p - q, m R )] 



(19) 



where S(p,m) = (p 1 + m)A(p, m) is the fermion propagator, A(p,m) being the free propagator for a scalar field 
of mass m and is given by 

A(p,m)= 2 ~\ ■ (20) 

Also included is a monopole form factor F(q) = A 2 /A 2 + q 2 with A = 2 GeV [12] to take into account the finite 
size of the pNR vertex. 

The corresponding expression for the case of loop graphs with spin 3/2 resonances is given by 



n^(g) =U F {^-) 2 j ■0^Tr[T^S(p,m N )T^S l3a (p- q ,m R ) 



(21) 



where the spin-3/2 propagator is S^ u (k,m) = (ljt + m)K^(k)(A(k,m) with K^ iu (k) = -g lxv + ^zk lx k v + \^ li ') v + 
3^(7/x^i/ — lukfj). Obtaining the vertex factors and T^ a from the interaction Lagrangians © and ([7]) both 
the expressions (|19|) and (|2ip can be expressed in the general form 

f d 4 p 

~R[iu{q) = iJ j—r^L flv {p,q)A(p,rn N )A{p-q,rn R ) (22) 

where the factor L^ iu (p,q) consists of the trace over Dirac matrices appearing in the two fermion propagators 
along with their associated tensor structures, isospin and form factors coming from the pNR vertex. Since the 
self-energy is transverse, L^ v can be expressed as 

lT(p, q) = I F 2 [a(p, q)A^ + /3(p, q)B^ + 7 (p, q)C^\ (23) 
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where the three gauge-invariant tensors A^ v , B^ v and C^ v are given by 

A nv{q) = -g^v + q^Qu/q 2 , 
Bftu(q,p) = q 2 p^Pu~ q ■ p(q»Pv + PtiQv) + {q-pfg^, 

C^(q,p) = q^P^Pu- q 2 {q-p){q^Pv+P l ,qv) + {q-pfq l ,qv ■ (24) 

The coefficient functions a(p,q), /3(p,q) and ^(p-.q) for the different loops are tabulated in the appendix. 

We now extend the vacuum self-energy to the nuclear medium. In the real-time version of thermal field 
theory that we are using the propagators assume the form of matrices. The spin and isospin structure of the 
self-energy graph remaining the same, it is only the scalar part A.(p, m) of the propagators that assumes a 
matrix structure. The required 11— component of the fermion propagator is given by 

E u {p) = A(p) + 2niN(p )5(p 2 - m 2 ); N(p ) = n + (uj)9(p ) + n„(u)9(-p ) . (25) 

The function n±(oj) = — - is the Fermi distribution where the ± sign in the subscript refers to baryons 

and anti-baryons respectively, uj = ^/p 2 + m 2 and fi is the baryonic chemical potential which is taken to be 
equal for all the baryons considered here. Expressed as 

nll 1 f 1 — n + n+ 1 — n_ n_ 



E 11 (p) = -—[ ^ + ^ ^—^ — (26) 

2uj \pq — uj + it p — uj — ie pq + uj — ie po + uj + ie J 

the first and the second terms can be identified with the propagation of baryons above the Fermi sea and holes 
in the Fermi sea respectively [30J while the third and fourth terms correspond to anti-baryons. 

As noted in the previous section, the in-medium self-energy function of the p can be obtained from the 
11-component of the thermal self-energy matrix. For the one- loop graphs shown in Fig. 1, the latter is given by 

/eft 
—^-L^{p, q)E n (p, m N )E n (p - q, m R ) . (27) 



Upon inserting the form of E 11 from (|25p we get three types of terms. One is the vacuum contribution 
involving the vacuum parts of the two propagators, the other two being medium dependent, one linear and the 
other quadratic in the thermal distribution function. Performing the po— integration and using the relations 
(|17p connecting the real and imaginary parts of the 11-component of the self-energy matrix with those of the 
diagonal element (defined using a bar), the self-energy function is written as 



^(qo,q) 



d 3 p 



(2tt) 3 4uj n ujr 



+ 



qo-uj N + UJR + ie(q )rj q + uj N -uj r + ie(q )r) 
Lf (l-n£)-L£V? Lf (l-n^)-Lfn* 



q - uj n - uj r + ie(q )r] q + uj N +uj r + ie(q )r) 



(28) 



where n N = u(ujn) with ujn = yP 2 + m%, n R = n(uj R ) with uj r = J (p — q) 2 + m 2 R and L^ u ,i = 1, ..4 denote 
the values of L^(po) for po = ujn, -Wjv, qo + ujr, qo — ujr respectively. 

Let us first consider the imaginary part of the self-energy. The retarded self-energy defined by (|18p can be 
easily read off from the self-energy function (|28l) to get 

d 3 p 



Imn-( go ,^ = -vrcoth(/W2) 7 (27r) 3 4a ^ 
[Lf {(1 - n£ - nX)5(q - uj n - uj r ) + (nj - n*)6(q -uj n + uj r )} 

+L^{{n R l - n N )5(q + uj n - ujr) - (1 - - n*)5{q Q + uj n + ujr)}} (29) 
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Irn(q„ 



Landau \/( m R — m^) 2 + q 2 



Unitary 



Inl R Z x 2 _l -2 Re ^ 

111 v /( m R + m N) 2 + q 2 

Figure 2: The positions of the Landau and unitary cuts in the complex plane for q 2 > 0. 



in which the factors Lg v 4 have converted to L^ u 2 respectively in association with the (5-functions. Following [161 [T7] 
it is interesting to relate the terms appearing in the above expression with scattering and decay processes 
involving the p, nucleon (iV) and the heavy resonances (R) in the medium. The delta functions in each of 
the terms in (|29p precisely define the kinematic domains where these processes can occur. The regions which 
are non-vanishing give rise to cuts in the self-energy function. Thus, the first and the fourth terms are non- 
vanishing for q 2 > (jur + m^) 2 giving rise to the unitary cut and second and third terms are non- vanishing 
for q 2 < (uir — tun) 2 giving rise to the Landau cut. Note that the unitary cut is present in vacuum but the 
Landau cut appears only in the medium. 

Consider, for example the first term with 1 — n+ — n R in fj29f) . Written as (1 — n+)(l — n?) — n+n R 
this indicates a process in which a (virtual) p decays into a NR~ 1 pair with the Pauli blocked probability 
(1 — tv±_){1 — n_) minus the process in which NR^ 1 pair gets absorbed in the medium with a statistical weight 
factor n+n R . This process can obviously take place for p's with invariant mass (yq 2 ) > (ttir + mjv)> a 
requirement that is in conformity with the kinematic threshold of the unitary cut coming from the associated 
5— function. 

The kinematic domains where the four terms contribute can be summarised as follows. For go > 0, the first 
term in (j28|) contributes for time-like values of q 2 , the second at space- like q 2 and the third at all q 2 . Likewise, 
for negative values of the variable qo, the second term is non-zero at all q 2 , the third at space- like q 2 and the 
fourth only for time- like q 2 . In view of the fact that the spectral function of the p will be measured in the 
invariant mass spectra of lepton pairs we will henceforth confine ourself to the kinematic region qo > and 
q 2 > 0. The position of the relevant cuts in the complex energy plane for this region are shown in Fig. [2] where 
we have ignored the portion of the Landau cut for qQ < |g|. (For a discussion on the branch cuts on the entire 
qo axis, see |17|). Furthermore, we will not include the unitary cut contribution in the imaginary part since 
the threshold of this cut begins at qo > tur + tun which being far away from the p pole is not expected to 
contribute to the p spectral function in a substantial way. Thus only the Landau cut as given by the third term 
in eq. ([29]) will be considered. 

Collecting the Landau contributions from both the diagrams using (jlip we finally write down the imaginary 
part of the p self-energy due to baryonic loops 



ImU^(q ,q) = vr J 



d 3 p 



(2tt) 3 2u n 



2u R > 
2ojr 



{(n^ - n+)6(q + uj n - ur>)} 

(30) 

where ujr/ = (p + q) 2 + m R- The two terms in this expression describe the contributions from scattering 

processes. The factor (rt+ — vS. ) expressed as (1 — n^')n^ — (1 — n^)n^' can be interpreted as the probability 
of a p meson scattering on a nucleon from the medium producing a resonance minus the process in which 
it scatters from the resonance to produce a nucleon, the final states in both cases being Pauli-blocked. The 
corresponding processes involving anti-baryons are included in the second term. 

Let us now proceed to evaluate the integral over the momenta in the NR~ l loop. The integral over cos# 

d 3 p = —2tt yu^f — m 2 N LdNduJNd{cos6) is done using the 5— function. Also the condition | cos 6*| < 1 puts 



m 
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restriction on the range of integration over cj^v- A substantial simplification is obtained by changing the 



integration variable to x using w/v 



-V 



( — Qo + \q\%) where S = q — m R + m N to get 



w 



dx L tlu (x){n + (q Q + uj N ) - n + (uj N ) + n_(g + uJn) ~ U-(uj n )} 



(31) 



IV 



where W = y 1 — 4g 2 m^ r /5' 4 and the factor L^ v in terms of the variable x has the same form for both diagrams 
in Fig. 1. This is shown in the appendix. 

The real part consists of principal value integrals which remain after removing the imaginary part from (|28p . 
Note that unlike the imaginary part, the real part of the self-energy at a given value of q receives contribution 
from all the four terms. 

Up to now we have been treating the baryon resonances R in the narrow width approximation. It is indeed 
necessary to consider the width of the unstable baryons in a realistic evaluation of the spectral function. For 
this, we follow the procedure (see e.g. [31|, I32j ) of convoluting the self energy calculated in the narrow width 
approximation with the spectral function of the baryons. This approach has the advantage that the analytic 
structure of the self energy discussed above remains undisturbed. 



1 r 



m R +2T R ^ 

dM-hn 

m R -2T R 7T 



M-m R + ^Tr(M) 



(32) 



with N R 



m R +2F R y 

dM-lm 

m R -2Y R 7T 



1 



and T R {M) = T R ^ N7r (M) + T R ^ Np (M). As a conse- 



M-m R + ±r R (M)_ 

quence of this convolution, the sharp ends of the regions of non-zero imaginary part smoothly go to zero at 
a higher value of M depending upon the width of the resonance. This is shown in Fig. [3] for the iV*(1520) 
resonance. 

Nuclear medium at finite temperature is also substantially populated by mesons which modify the p propa- 
gation in the medium in a non-trivial way. This has been studied [T7] following the same procedure as described 
here for mesonic loop graphs with one internal pion line and another meson line h where h = 7r, cj, hi, a\ using 
interactions from chiral perturbation theory. Collecting the real and imaginary parts, the self-energy from 
mesonic loops can be written as 



n^/(<?o,£) = 
Nf(l + n 



d 3 k 



1 



N^n h 



qo - uj-k + + ie(qo)v qo + - Uh + ^(<?o)?7 



(27r) 3 4u) n Uh 



+ 



q +uj n -uj h + ie(q )t] q + uj n + co h + ie(q )r] 



(33) 



where the Bose distribution functions n n = n(a; w ) with lj„ = y k 2 + m 2 and nh = n(ujh) with Uh = y (q — k) 2 + m^. 
N^ u {i = 1,4) are the values of N fJ,u (ko) for ko = co n , —uj n , qo + ujh, qo — respectively and can be expressed in 
terms of the gauge- invariant tensors (|24h . The complete expressions are provided in the appendix. Using the 
procedure described above we have improved upon the calculations in [17] by including the width of the heavy 
mesons a\ and h\. In this case we use a slightly different formula 



W^(q;m h ) 



1 



with Nh 



(m h +2T h ) 2 i 

dM 2 -Im 

{m h -2T h y 7T 



{m h +2T h y 1 

dM 2 -Im 

{m h -2T h f K 
1 



1 



M 2 -m 2 + iMT h {M) 



(34) 



and T h (M)=T h ^ p7T {M). 



_M 2 -m 2 + iMT h (M)_ 

The transverse and longitudinal components can then be obtained from the self-energy tensors using the 
relations (fTBT). 
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Figure 3: Imaginary part of p meson self-energy showing the individual contributions for different NR loops. 
Left panel shows results for q = and the right panel shows the transverse (solid) and longitudinal (dotted) 
parts for q = 300 MeV. The long dashed line in the upper panel on the left shows the imaginary part of the 
iViV* (1520) loop evaluated in the narrow width approximation. 
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Figure 4: Same as Fig. [3] for the real part 



5 Numerical Results 



In this section we present the results of numerical evaluation beginning with the imaginary part of the p self- 
energy as a function of the invariant mass \fq* = M for two values of the three momentum. The Landau 
cut contribution starts from M = for all values of the three momentum. Shown in Fig. [3] left panel are 
the contributions from the individual NR loops for a p meson at rest. The A r iV*(1520) loop makes the most 
significant contribution followed by the iV*(1720) and A(1700). The right panel shows the corresponding results 
for q = 300 MeV where the transverse and longitudinal components IT and q 2 Hi have been shown separately 
by solid and dotted lines respectively. (Note that for a p meson at rest IT = <7q IT; . ) The corresponding results 
for the thermal contribution to the real part are shown in Fig. HJ The divergent vacuum contribution in this 
case is assumed to renormalize the p mass to its physical value. Also shown by the long dashed lines in the left 
upper panels of Figs. [3] and H] are the corresponding contributions to the real and imaginary parts coming from 
the AW* (1520) loop computed in the narrow width approximation. 
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M (GeV) M (GeV) 

Figure 5: (Left) Imaginary and real parts of p meson self-energy showing the individual contributions for 
different it — h loops. (Right) The total contribution from meson and baryon loops. 




M (GeV) M (GeV) 

Figure 6: The spectral function of the p meson for (left) different values of the temperature T and (right) 
different values of the three-momentum q. 

Next we show the results of the spin-averaged p self-energy defined by 

n = i(2n t + q 2 n{) . (35) 

where the transverse and longitudinal components are obtained using (I16p . In Fig. [5] left panel the imaginary 
and real parts of tt — h loop graphs are shown in the upper and lower panels respectively. The Landau and 
unitary cut contributions for the it — uj loop are clearly discernible though the contribution at the p pole is 
dominated by the a\. On the right panel we plot the total contribution from the baryon and meson loops for 
two values of the baryonic chemical potential. The small positive contribution from the baryon loops to the 
real part is partly compensated by the negative contributions from the meson loops. The substantial baryon 
contribution at vanishing baryonic chemical potential reflects the importance of anti-baryons. 
We now turn to the spin averaged spectral function given by 

ImG( ? ) = \{2lm.G t + qHmG^ (36) 
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Figure 7: The spectral function of the p meson for different values of the baryonic chemical potential p 
where 

Im G t i (q) = _ tnt U -r-f— . (37) 

{M 2 _ m 2_ (1 ^ 2) £ Re 11^)2 + {(l,5^knjp 

First, in Figj6]left panel we plot the spectral function at fixed values of the baryonic chemical potential and three- 
momentum for various representative values of the temperature. We observe an increase of spectral strength 
at lower invariant masses resulting in broadening of the spectral function with increase in temperature. This is 
purely a Landau cut contribution from the baryonic loops arising from the scattering of the p from baryons in 
the medium. However, we do not observe much variation with the three-momentum of the p as seen from the 
figure on the right panel. 

We then plot in Fig. El the spectral function for various values of the baryonic chemical potential for a fixed 
temperature. For high values of p we observe an almost flattened spectral density of the p. 

Owing to differences in the various approaches to the evaluation of the p spectral function followed in the 
literature, a direct numerical comparison with earlier results does not appear to be meaningful. These differences 
are at the level of the basic formulae arising from the type of couplings of the nucleon and the rho fields with 
the various baryon resonances considered as well as in the form of the propagators used in the calculations. 
There also exist differences at the level of formalism employed in the evaluation of the p self-energy. We thus 
end this section by showing how the in-medium spectral function of the p is manifested in the dilepton emission 
rate. This rate from thermalised hadronic matter is given by |33j 

= ~3^ H{M2)nBE{qo) 9^T^(q , q) (38) 

where ImT^ is the imaginary part of the (retarded) two-point function of vector currents which can be obtained 
from © using relations analogous to flU])- The quantity H(M 2 ) = (1 + 2mf/M 2 ) (1 - 4ra 2 /M 2 ) 1 / 2 is of the 
order of unity for electrons and will be omitted henceforth. In the low invariant mass (M) region, ImT^ is 
usually expressed as a sum over the spectral densities of the vector mesons p, ui and 4> pp. This is however 
justified only in vacuum. The vector mesons p and uj can in general undergo mixing in the presence of matter 
which can lead to non-trivial modifications, for example, of the electromagnetic form factor of the pion [3j| and 
consequently the invariant mass distribution of lepton pairs. In the following, we will consider only the p pole 
contribution which is known to play the most dominant role. The dilepton rate in this case can be expressed 
in terms of the spin-averaged spectral function of the p (|36p getting, 

^j^n BE (q )K p lmG(q ,q) (39) 



dM 2 qT<lqTdy tt 2 
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M (GcV) 

Figure 8: The lepton pair emission rate at T = 175 MeV with and without baryon (B) loops in addition to the 
meson (M) loops 

where we have used = K p G^ u as defined earlier. Integrating over the transverse momentum and rapidity 
y of the electron pairs we plot dR/dM 2 vs M in Fig. [8] for T=175 MeV. Because of the kinematical factors 
multiplying the p spectral function the broadening appears magnified in the dilepton emission rate. A significant 
enhancement is seen in the low mass lepton production rate due to baryonic loops over and above the mesonic 
ones shown by the dot-dashed line. The substantial contribution from baryonic loops even for vanishing chemical 
potential points to the important role played by antibaryons in thermal equilibrium in systems created at RHIC 
and LHC energies. 

6 Summary and Discussion 

We evaluate the p self-energy to one loop in nuclear matter at finite temperature and baryon density. Loop 
graphs involving the nucleon and 4-star iV* and A resonances up to spin 3/2 were calculated using gauge 
invariant interactions in the framework of real time thermal field theory to obtain the correct relativistic 
expressions for the p self-energy. The singularities in the complex energy plane were analysed and the imaginary 
part obtained from the Landau cut contribution. Results for the real and imaginary parts at non-zero three- 
momenta for various values of temperature and baryonic chemical potential were shown for the individual loop 
graphs. Adding the contributions from mesons obtained in the same formalism, the spectral function of the p 
was observed to undergo a significant modification at and below the nominal rho mass which was seen to bring 
about a large enhancement of lepton pair yield in this region. 

It may be emphasised that the determination of the p spectral density at finite temperature and baryon 
density by an explicit evaluation of loop graphs using thermal field theoretic techniques such as performed here 
is of relevance in view of precision data from experiments at RHIC and LHC as well as from the FAIR facility 
at GSI in future. But an actual comparison with data will involve a space-time evolution of the static rates 
using a framework like relativistic hydrodynamics. Efforts in this direction are in progress and will be reported 
in due coarse. 
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Figure 9: Contour in the complex time plane for real time formalism 



8 Appendix 

8.1 The real-time propagators and the self-energy matrix 

In the real time formulation of thermal field theory a two-point function of local operators assumes a 2 x 2 
matrix structure on account of the shape of the contour in the complex time plane. Thus the current correlator 
T^ v [35] as well as propagators for various fields [20J assume a matrix structure. For the contour shown in Fig. [9] 
the components of a free scalar propagator are as 



D 



ii 



{D 



22\* 



A(q) + 2mn5(q 2 - m 2 ) 

2\ 



D 12 = D 21 = 2TTiy/n(l + n)5(q 2 
where A(q) is the Feynman propagator in vacuum, 



m 



q 2 — m 2 + ie 



The thermal propagator may be diagonalised in the form 



D a \qo,q) = U ac (q )[di a g{A(q ,q),-A*(q ,q)}rU ab (q ) 



cd t rdb 1 



(A.l) 
(A.2) 
(A.3) 



with the elements of the diagonalising matrix as 

U u = U 22 = y/l + n, 
Using the (transverse) vector propagator given by 



U 



12 



u 



21 



n . 



(A.4) 



in the Dyson equation and the fact that U diagonalises not only the free propagator, but also the complete 
one \20\ I21| it turns out that the self-energy matrix IP?* is also diagonalisable by (U~ 1 ) ab , 



n#(?) = [c/- 1 (9o)] ac [diag{n^(g),-n^( g )}] C(i [c/- 1 (9o)] d6 • 



The relations 



Re n 
ImTT 



[IV 



Ren 



ii 



[IV 



e(g )tanh(/3g /2)Imn 



ii 



(A.5) 



(A.6) 



trivially follow showing that the thermal matrices are actually given by a single analytic function which essen- 
tially coincides with the corresponding result from the imaginary time formulation. 
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R 


J p 


f 


OL 





1 


iV(940) 


1 + 

2 


7.7 


\ a l/2+ 


2&1/2 


hi/2 


iV*(1520) 


3~ 
2 


7.0 


"3/2- 


03/2 


73/2 


iV*(1650) 


1- 
2 


0.9 


"1/2- 


01/2 


7l/2 


iV*(1720) 


3 + 
2 


7.0 


«3/2+ 


03/2 


73/2 


A(1232) 


3 + 
2 


10.5 


«3/2+ 


03/2 


73/2 


A(1620) 


1- 
2 


2.7 


"1/2- 


01/2 


71/2 


A (1700) 


3- 
2 


5.0 


"3/2- 


03/2 


73/2 



Table 1: Table showing the coefficients a, and 7 for loops containing the various resonances considered. 

8.2 The factor L^ u in baryonic loops 

The factor L^ u which appears in the loop integrals is given by 

L^(p,q) = I F (HM] [aA^ + 0B^ + jCn (A.7) 
\ m P J 

where A^ u , B^ u and C^ v are the transverse tensors defined in (|24p . The coefficient factors a, 0, 7 are given in 
Table 1 where 

"1/2+ = 4g 2 (p 2 - m N m R + p- q) 
"1/2- = a 1/2+ (mjv -> -mjv) 

A/2 = 8 

7i/2 = 

"3/2+ = J^2"[P 2 (P 2 - 3^ 2 ) +p- q(3p 2 + q 2 ) 

+3m 2 R (p 2 + 2m N mn + p ■ q) 
-2m N m R {p 2 + q 2 - 2p ■ q)] 
"3/2- = Oi 3 / 2+ (m N -> -mjv) 

03/2 = 4(1+^/3771^) 

73/2 = "4/3m|. (A.8) 

The corresponding expressions of a, and 7 for the second graph of Fig. 1 can be obtained by replacing 
q — > —q as indicated in (jlip . The factor 1/2 in front of the coefficients for the NN loop is put to prevent a 
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double counting of this contribution. The coupling constants / have been obtained as in \13\ I15j. The values 
are also in reasonable agreement as seen in Table 1. 

A considerable simplification in the imaginary part can be achieved in the sum of the two diagrams in Fig. 1 
by a change of variables. Note that only the factor p ■ q which appears with various powers in the tensors 
A 11 ", B^ v and C^ v and in the factors a takes on different values in the various terms in the expression for the 
imaginary part of self energy. We recall for convenience the imaginary part coming from first diagram in Fig. 1, 

(i), f d 3 p 



lmU {1 \q ,q) = -7r J 



(2tt) 3 4:U; n u}r 

[Li(q){(l - - n*)5(q - u N - ujr) + (nj - n+)%o - oj n + ujr)} 

+L 2 (q){(n B : - n^)% + - ur) - (1 - r£ - n*)5(q +uj n + cjr)}] (A.9) 

where we have dropped the Lorentz indices for brevity. For the first two terms po = ojn and on integration over 
the angle using either of the two delta functions one has p-q = — 1(5* 2 — 2qouiN) with S 2 = q 2 — m 2 R + m 2 N . One 
thus gets p ■ q = poqo — p '■ q = Changing variable to x using ojn = j^(Qo + \o\ x ) one nas 

A» = -3, A 00 = ] \ (A.10) 



B» = mlq 2 + —, Boo = -^(l-x 2 ) (A.ll) 
C£ = q\mW ~ C 00 = ^x 2 . (A.12) 



For the last two terms, p ■ q = — \{S 2 + 2go^7v) and po = —ojn so that here too p ■ q = 4^. Defining x in this 

2? 



case through ui^ = -§^(—qo + \o\ x ) one recovers the same expression for the tensor component as in the first 



two terms. 

The corresponding expression for the imaginary part from the second graph of Fig. 1 is obtained by changing 
q — > — q in the expression for Imll given by (|A.9jl getting 

Im5(2 W-) = -7r/ ' p^ * 

[L 2 (-q){(l -n^ - n+)6(q - u) N - u R >) + (n^ - n^')S(q - ui N + u R >)} 

+L 1 (-q){(n* - n")5(qo + w N - u a ) - (1 - n? - n R ')8(q + uj N + wr>)}] (A.13) 



where ojr> = y (p + q) 2 + Tn R ,. In this case, for the first two terms p ' ■ q = \{S 2 — 2qoU>N) which along with 

Po = —ujn gives p ■ q = —^2~- The same is obtained for the last two terms for which p '■ q = \{S 2 + 2qoU}^) and 
Po = ujn- Making a change of variables as before, identical values of the tensor components are obtained as in 
(|A.12|) as a consequence of the fact that the gauge invariant tensors A^ u , B^ v and C M „ are even under q — > —q. 
The coefficients a, (3 and 7 also remain unchanged in the two diagrams the sign of p ■ q in the expressions 
remaining the same under the combined effect of q — > —q and a reversal in its mag nitude (£- -)■ -£-). The 
value of L{x) thus comes out to be the same for all cases and we end up with the final expression given by (|31|h 
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8.3 Expressions for N^ u for mesonic loops 

The expressions for N^ u appearing in the p self-energy (|33|) for ir — h loops have been obtained in [T7j and are 
given below, 



Nj£\q,k) = -4(|r) (iV + g 2 ^, 



iVj^fak) = - 2 (|-) 2 (^--i-^) (A.14) 

where the constants G p = 69 MeV, F n = 93 MeV, g\ = 0.87, 52 = 1-0 and 53 = 1.1. These can be simplified 
as shown above and finally expressed in terms of the tensor components defined in (|A.12j) . As in the case of 
baryon loops the factors N^u have also been multiplied by the square of the monopole form factor F(q) defined 
earlier. 
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